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EVALUATION 


The  necessity  for  uore  complex  software  systems  In  such  areas  as 
command  and  control  and  avionics  has  led  to  the  desire  for  better 
methods  for  predicting  software  errors  to  insure  that  software 
produced  is  of  higher  quality  and  of  lower  cost.  This  desire  has  been 
expressed  in  numerous  Industry  and  Government  sponsored  conferences, 
as  well  as  in  documents  such  as  the  Joint  Commanders'  Software 
Keliabillty  Working  Group  Report  (Nov  1975).  As  a result,  nunerous 
efforts  have  been  initiated  to  develop  methods  for  determining  the 
optimal  policy  for  maintenance  of  an  operational  software  system, 
however,  early  efforts  have  not  developed  any  consistent  or  generally 
applicable  software  maintenance  policy. 

This  effort  was  Initiated  in  response  to  this  need  for  developing 
better  and  more  accurate  software  error  prediction  models  and  fits 
into  the  goals  of  PADG  TPO  No.  5,  Software  Cost  Reduction  (formerly 
PAf)C  TI’O  Vo.  II,  Software  Sciences  Technology),  in  the  subthrust  of 
Software  duality  (Software  Voile  1 tng) . This  report  summarizes  the 
development  of  a haveslan  methodology  for  determining  the  optimal 
policy  for  maintaining  an  operational  software  system.  The  importance 
of  this  development  is  that  it  represents  the  first  attempt  to  develop 
operational  software  maintenance  policies  that  more  closely  reflect 
llit-  actual  software  error  detection  and  correction  process. 

The  theory  and  equations  developed  under  this  effort  will  lead  to  much 
needed  predictive  measures  for  use  by  software  maintenance  personnel 
in  providing,  better  and  more  efficient  maintenance  of  operational 
software.  In  addition,  the  associated  confidence  limits  and  other 
re  la  ted  statistical  quantities  developed  under  this  effort  will  insure 
more  widespread  use  of  these  modeling  techniques.  Finally,  the 
predictive  measures  and  equations  developed  under  this  effort  will  be 
applicable  to  current  Air  Force  software  development  projects  and  thus 
help  to  produce  the  high  quality,  low  cost  software  needed  for  today's 
systems . 


ALAN  N.  SUKKRT 
Project  Engineer 


1.  INTRODUCTION 


\ 


In  this  report  we  discuss  th"  problem  of  determining  an 
optimum  correction  limit  policy  for  a large  software  system  which 
is  subject  to  random  occurrences  of  errors.  When  in  error  occurs, 
a corrective  action  is  undertaken  to  remove  it.  Such  an  action 
can  be  scheduled  at  two  levels,  which  we  call  Phase  I and  Phase  II. 
By  Phase  I we  mean  that  the  corrective  action  will  be  undertaken 
by  the  programmer  while  Phase  II  action  is  undertaken  by  a system 
analyst  or  system  designer.  First,  Phase  I corrective  action  is 
scheduled  for  a specified  time  T.  If  the  error  is  not  corrected 
in  this  time,  it  is  referred  to  Phase  II.  This  sequence  of  cor- 
rective actions  in  an  operational  phase  is  shown  in  Figure  1.1. 

Our  objective  is  to  determine  the  optimum  value  T*  of  T which 
minimizes  the  long  run  average  cost.  Two  models  are  developed 
for  this  purpose.  In  the  first  model  (Section  2)  we  assume  that 
the  cost  cf  observations  of  error  occurrence  and  correction  time, 
prior  to  the  implementation  of  the  optimum  policy,  is  negligible, 
"'he  second  model  (Section  3)  incorporates  the  cost  of  observations. 


2.  MODEL  WHEN  COST  OF  OBSERVATIONS  IS  NEGLIGIBLE 


The  following  assumptions  are  made  for  model  development: 

(i)  The  error  occurrence  time  in  a software  system  has  an 
exponential  distrib_  .ion  with  an  unknown  mean  \ . 

(ii)  The  error  correction  time  at  Phase  I is  exoonential  with 
an  unknown  mean  ^ . 

(iii)  The  Phase  II  error  correction  time  has  a general  distri- 
bution with  a known  mean  ^ • 

(iv)  Appropriate  prior  distributions  can  be  chosen  for  \ and 
“1  * 

2.1  Predictive  Distributions  of  Y and  X 

Let  random  variables  X and  Y denote  the  error  occurrence 
time  and  Phase  i enor  correction  time,  respectively.  From  assump- 
tions (i)  and  (ii)  the  probability  density  functions  are 

f(x|\)=-£eX,X  x>0,  \>0  (2.1) 

g(yiu,)  = -7-  e 1 y>0,  ^>0  (2.2) 

In  this  report,  we  develop  suitable  expressions  by  consider- 
ing the  conjugate  priors  for  ^ and  X which  are  inverted  gamma 
distributions  given  by 

P(nL) 

ana 

P(X) 
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r(ftl)  *1 


“ ( Of  j + 1 ) ~ ® 1 ' ^ 1 


a 1 ' ® 1 > 0 ' 


(2.3) 


e2  -(a2+i)  *e2/x 

n» ,)  x 


o2 , b2  > 0 


(2.4) 


The  expressions  for  any  other  reasonable  priors  for  and 

X can  be  developed  similarly. 

Also,  let  x = (x1,x2,  . . . ,xn)  and  X.~^yl'y2 yn)  be  the 

observed  values  of  n error  occurrence  times  and  n error  correction 
times,  respectively. 

Now  we  obtain  expressions  for  the  predictive  distributions 

of  V and  X and  also  obtain  the  Bayesian  estimates  y , and  x , 

n+  1 n+ 1 

which  will  be  used  to  obtain  the  cost  function. 

For  given  observations  ^ , the  likelihood  function  of  is 


y>i 

e 1-1 


(2.5) 


The  posterior  distribution  of  ^ is  obtained  from  Bayes 


theorem: 


P(^ly)  = 


1 !y)p^i) 

J )&l 


(2.6) 


Substituting  the  expressions  for  p(^)and  f(^ly)  from  (2.3)  and 
(2.5),  we  get  the  posterior  distribution  of  ^ as 


p(ui'X)  * 


_ (Sl  i=i  yi>  -(0-,+n+l)  ~*(B1  + “ yi)/u1)  (2.7 ) 

V — * >1,  e 1 1 

r(wj_  + n) 

Using  this  posterior,  the  predictive  distribution  g(yly)  of 


a^+n 


error  correction  time  at  Phase  I is  given  by 


g(yly>  - JQ  9(yU1)p(^1lx)d^i 


(2.8) 
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-V 


Substituting  the  expressions  for  g(ylp^)  and  p(^ly)  from  (2.2) 
and  (2.6),  respectively,  we  get 


«<*,*>  - 


- (n+cr^+l) 


(2.9) 


E y • + 0,  E y . + P. 

i=l  i=l  1 1 


The  cumulative  predictive  distribution  to  some  specified  time  t is 


G(t|y)  = J g(yly)dy 
~ 0 ~ 


“ 1 " (~n  " ) 


~ (n+o-^) 


(2.10) 


i=l  1 x 


We  define  the  predictive  phase  I error  correction  rate  as 


g(t ly) 

r(t|  y)  = -- 

G(tl^) 


(2.11) 


so  that 


(2.12) 


s yi  + Pi 

i=l  1 A 


2 Yi+Pl 

i=l  1 


where  G (t)  sl-G(t)  . 


The  predictive  distribution  f(x|x)  of  the  error  occurrence 
time  can  be  similarly  obtained. 


2.2  Bayesian  estimates  x_ and  y . 

— * n+ 1 n+ 1 

From  the  predictive  distributions  of  X and  Y the  Bayesian 

estimates  x . of  the  time  to  (r.+  l)st  error  occurrence,  for 
n+ 1 

given  x and  the  (n+l)st  error  correction  time  for  given  y are 
easily  obtained,  since 


5 


Similarly 


where 


J*  G(tl^)dt 

* yi  + Bl 


oL  + n-1 


, = f F(tlx)dt  , 


F (tl  x)  = f f (x  |x)dx 
~ J 0 ~ 


(2.13) 


(2.14) 


(2.15) 


which  is  a cumulative  predictive  distribution  to  some  specified 


time  t . Hence 


£ *.  + e2 

i=l  1 z 
a 2 + n-1 


(2.16) 


2 . 3 Cost  Function 

Let  c^  (c2)  be  the  cost  per  unit  time  of  error  correction  in 
Phase  I (Phase  II)  and  the  costs  be  linear  functions  of  time.  From 
assumption  (iii)  Phase  II  error  correction  time  has  some  arbitrary 
general  distribution  with  a known  mean  If  we  consider  one  cycle 

to  be  the  time  from  the  beginning  of  (n+l)st  operation  to  the 
beginning  of  (n+2)nd  operation,  then  the  expected  cost  in  one  cycle 


E(C)  » c1  J^G(tl^)dt+  c2vi2G(Tly)  , 


(2.17) 


where  T denotes  the  scheduled  correction  limit  time  in  Phase  I. 
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The  expected  length  of  one  cycle  is 


E (h)  - xn+1  J*  0(tly)dt  + n2  G(Tly)  , 

and  hence  the  long  run  expected  cost  per  unit  time  iu 

E(C) 

C(T)  ' "eTlT 


(2. IB) 


C (T)  - 


cl  S Q(tly)dt  + c2  vl2  Q(Tly) 
__  _ 

o(tl£)dt  ♦ p2  b(TlY) 


(2.19) 


This  is  the  cost  function  which  we  want  to  optimize  to  obtain  the 
optimum  policy. 


From  (2.19),  we  note  that 


C(0)  * 7 


C2  ^2 


(2.20) 


C(«)  - 


(2.21) 


Xn+ 1 * Vn+  1 


where  f . is  the  Bayesian  estimate  of  y for  given  data  y . 

ITT  X 

Also,  note  that  T - 0 means  that  the  errors  are  corrected  only 
at  Phase  II  while  T*«  means  tha*  they  are  corre  ted  at  Phase  I. 

To  obtain  an  optimum  T*  which  minimizes  the  long  run  average 
coat  per  unit  time,  C (T)  , we  need  the  following  theorems  and 
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corollary.  Theorems  2.1  and  2.2  are  the  special  cases  of  the 
theorems  proved  in  Appendices  A and  B respectively. 

Theorem  2 . 1 

Assume  c^  < c2  . Then  under  the  following  condition 


r(Oly)  > 


cl(xn+l  + tt2)  ~ c2  ^2 
c2  xn+l  ^2 


(2.22) 


there  exists  a finite  and  unique  T*  which  satisfies 


r(Tly)(c2  xn+1+  (Cj-c^  J*  G(tly)dt) 


+ (Cj-c^GtTly) 


°1  xn+l 


(2.23) 


Theorem  2 . 2 


If  the  above  conditions  are  satisfied  then  there  also  exists 


a finite  and  unique  upper  bound  T(>T*)  such  that 


r(Tly) 


ClXn+l 


l*2(c2Xn+l+  (C2_C1)  yn+lJ 


(2.24) 


This  upper  bound  can  be  used  to  obtain  an  initial  value  for 
solving  the  nonlinear  equations  in  T*  . 

Corollary  2.1 

If  there  exists  an  optimum  T*  , then  the  associated  cost 
function  is  given  by 


c.-c-  r (T*  ly) 
C(T*>  * 1 ~ ^2  r (T*  iy) 


(2.25) 
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2. 5 Numerical  Example 


We  use  simulated  data  in  this  example  to  illustrate  the  cal- 
culations and  nature  of  various  quantities  in  the  determination 
of  T*  . 

Let 


cx  - 8000  c2  * 9000 

0^=0  02  = 0 

o2  * 0 P2  ■=  0 

H2  ■ 0.7 

The  simulated  data  (xn#yn)  are  9iven  i-n  Table  2.1.  Suppose 
n»10  data  points  are  available.  The  Bayesian  estimates  of  and 

y^^  are  obtained  from  (2.16)  and  (2.13)  as  X22-59.6O  and  yj^0*?8* 
respectively.  Such  values  for  various  n are  given  in  Table  2.2. 
For  the  case  n * 10  we  see  that  the  optimum  correction  limit  time 
is  T*  * 0. 90  hours  and  the  corresponding  minimum  cost  rate  is 
C (T*)  *99.44  dollars  /hour. 

Thus,  for  this  set  of  data,  we  will  schedule  corrective 
action  in  Phase  I for  0.90  hours  and  if  it  cannot  be  completed 
in  this  time,  the  software  system  will  be  referred  to  the  system 
analyst  for  corrective  action. 


TABLE  2.1 

Simulated  Values  of  x and  v 

n n 


TABLE  2.2 

Calculation  for  the  Optimum  Correction  Time  Policy 


xi»l(hr*) 


89.17 

121.74 

86.02 

67.23 

60.85 

63.70 


59.60 


Vl(hr*> 


12 

53.88 

13 

52.00 

14 

55.46 

15 

57.09 

16 

61.58 

17 

57.76 

18 

55.09 

19 

56.78 

20 

55.03 

2.98 

1.92 

1.36 

1.02 

0.88 

0.80 

0.90 

0.84 

.78 

.80 

.82 

.80 


T*  (hr.) 


C (T* ) 


0.32 

101.05 

0.73 

80.11 

0.02 

90.78 

0.38 

100.60 

0.90 

99.44 

0.66 

102.87 

0.50 

113.82 

0.68 

116.85 

1.45 

103.80 

0.75 

106.51 

1.02 

97.47 

1.45 

101.22 

2.16 

101. 14 

109.61 

112.97 


2.6  Sensitivity  Analysis  of  the  Optimum  Correction  Time  T* 


To  study  the  sensitivity  of  T*  to  changes  in  various  para- 
meters, we  look  at  the  plot  of  T*  versus  average  correction  time, 
y - E y^n  . Plots  for  various  values  of  c^,  c2>  a^,  c»2, 

02,  \,  ^ and  ^ are  given  in  Figures  2.1  to  2.4.  in  Figure  2.1, 

is  varied  while  other  parameters  are  kept  constant.  The  effect 
of  changing  c2  while  keeping  other  factors  constant  is  shown  in 
Figure  2.2.  Effects  of  changing  u2  ancl  n are  shown  in  Figures  2.3 
and  2.4,  respectively.  The  following  observations  can  be  made  from 
these  figures. 

(i)  T*  increases  with  a^  for  fixed  y and  the  slope  of  T* 
vs  y lines  is  independent  of  (Figure  2.1) 

(ii)  T*  increases  with  c2  for  fixed  y and  the  slope  of  T* 
vs  y line  is  independent  of  c2  (Figure  2.2) 

(iii)  T*  increases  with  f°r  fi-xecl  y and  the  slope  of  T* 
vs  y line  is  independent  of  (Figure  2.3) 

(iv)  T*  increases  with  n for  fixed  y but  the  rate  of  increase 
decreases  with  y (Figure  2.4). 


AVERAGE  CORRECTION  TIME 


OPTIMUM  LIMIT  TIME 


AVERAGE  CORRECTION  TME 


3.  MODEL  FOR  NONZERO  COST  OF  OBSERVATIONS 


In  this  section  we  develop  the  cost  model  for  the  optimum 
correction  limit  policy  by  incorporating  sampling  cost.  It  is 
assumed  that  the  sampling  cost  is,  a linear  function  of  the  sample 


Let  c be  the  sampling  cost* at  each  state  and  Cn(Tn)  be  the 
expected  cost  per  unit  time  unAl  the  completion  of  (n+l)st  correc- 
tive action  under  the  limit  cim»i  Tn  . Having  taken  n observations, 
if  we  decide  to  take  another  observation,  i.e.,  the  (n+l)st  observa- 
tion, then  C . (T  , ) is  the  cost  per  unit  time  until  the  completion 
n+ 1 n+ 1 

of  the  (n+2)nd  corrective  action  under  the  limit  time  Tfi+^  . 

3.1  cost  Function 

Let  the  length  of  the  nth  cycle  be  the  time  from  the  beginning 
of  nth  operation  to  the  end  of  nth  corrective  action.  Then  the 
expected  cost  at  the  end  of  (n+l)st  cycle,  given  n observations. 


E(Cn)  = nc  + c1  J*  G (t|  y)dt  + c2  G (Tnl  £) 


(3.1) 


The  expected  time  to  the  end  of  the  (n+l)st  cycle  is 


E|Ln>  ‘ iJi  *i  + *n+l*if1  *W0  0(tlJj)dt+»20(Tnlz) 


(3.2) 


The  expected  cost  per  unit  time  at  the  end  of  (n+l)st  cycle  is 

•« 

then  given  by 
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C„(T  ) 


E(C  ) 


n n E(L 


cn(V 


nc  + c1  G(tl  ^)dt  + c2  n2  G(Tnl 

i=1  Xi  + *n+l+  yi  + ;oD  °(tli:)dt+^25(Tnly) 


(3.4) 


I f we  decide  to  consider  the  next  cycle,  then  the  cost  rate 
function  to  the  end  of  (n+2)nd  cycle  is  similarly  obtained  as 


(n+l)c+clv*  G (tl  y)dt  + c.  .2  G (Tn+1l 


* Xi+  2xn+ 1 + ,r,  *i+ *1*1  G(t(>:)dt+  ^Gd-^iH;) 

1=1  1=1  0 


3.2  Optimum  Policy 

Our  objective  is  to  determine  the  optimum  sample  size  n* 
and  the  optimum  correction  limit  T*  such  that 


cnIT;’  s Cd.lITS.l>  • 


(3.6 ) 


Given  that  n observations  have  been  taken,  the  following  steps 
summarize  the  procedure  of  determining  these  quantities: 

(i)  Calculate  Cn(T£)  and  Cn+1(T*+1) 

(ii)  If  C (T* ) £ C_ , (T*  , ) , then  stop  taking  observations 
n n n+  i n+ 1 

and  employ  n*  and  T*  as  the  optimum  policy. 

(iii)  If  C (T*)  > C , (T*)  , take  (n+l)st  observation,  i.e. 
n n n+ 1 n 

let  n*n+l  and  go  to  step  (i). 


(3.5) 
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The  following  theorems  are  used  in  determining  c (T* ) and 

n n 

Cn+^(T*+^)  . Theorems  3.1  and  3.2  are  the  special  cases  of  the 
theorems  in  Appendices  A and  B,  respectively. 

Theorem  3 .1 

Suppose  c^  < c^  and  A a 0 . Then  there  exists  a unique  and 
finite  T*  satisfying 

Tn  _ 

r (Tniy)  (A+  (02-Cj^)  g (t  ly)dt)  + (c^c^G  (TR  ly)  « B (3.7) 

where 


A = c 


nc 


{i-iXl**n+l+wyi}  ‘ 

B‘^[el{islxi+*»l+isl»i}-"c]  • 

Also,  the  associated  cost  rate  function  is  given  by 

Cl-C2^2r(Tn) 


(3.8) 


(3.9) 


Cn<T£> 


1-‘*2r(Tn)  * 


(3.10) 


Theorem  3.  2 

If  the  conditions  of  Theorem  3.  1 hold,  then  there  exists  a 
finite  upper  bound  Tn  (>T*)  such  that 


r<V*>  • A+  <^>W 


(3.11) 
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1 


Th«Qr«n>  J.  J 

If  both  T*  ami  T*+  ^ exist,  then  the  following  relation- 
ship ho Ida t 


VTS' 

*>  c"‘l<T 

<*•> 

r(T*) 

,s,  r,T*»>' 

<c«> 

T*  > 

T*  , 

h (S) 

n+  1 

Ke lat lonah tps 

given 

in  Theorem 

( « . i : ) 


Appendix  C. 

t . J Numerical  Kx Ample 

In  this  example  we  use  emulated  data  to  illustrate  the 

determination  ot  n*  and  T*  . Simulated  valuers  ot  n and  v 

n n n 

for  vai  icus  n air  given  in  Table  l.l.  Suppose  the  values  or 
vai ious  guant it  lea  ale  as  fcllowSt 


C1  ' 

*1  * 

a , ” 

i*-»  * 


HOOD 

.8 

0 

0.7 


o , 

f ( e l 

f , * 0 


>1000 


•Id 


Then  the  values  of  *n+\>  ^1*  Tn*  Cn(Tn)  * Tnt  1 A,ul 

, ,T*  \ are  obtained  from  the  above  expressions  and  are  > nven 

v n*  l ' n+  1 

in  Table  J.  J.  From  this  table  we  see  that  for  n-  11, 

C'u  (T^>  - 21.74  and  (T^j)  * ■*'•*>*  so  that  0 1 1 ^ ' 1 l ' * •' ' 

Therefore,  the  optimum  policy  is  n*  - 1 1 aiul  T*-0.0'». 


JO 


TABLE  3.1 

Simulated  Values  of  x and  y 

n Jn 


n 

Vhr) 

*n<hr) 

n 

x (hr) 
n 

Vn0'r) 

1 

32.25 

1.69 

10 

3.72 

0.  15 

2 

34.77 

0. 12 

11 

50.85 

0.07 

3 

63.92 

0.23 

12 

64.89 

0.12 

4 

21.03 

0.41 

13 

0.76 

0.  29 

5 

39.42 

0.20 

14 

87.45 

1.33 

6 

9.91 

0.  37 

15 

64.12 

0.77 

7 

3.69 

0.  22 

16 

30.98 

1.37 

a 

2.42 

1.75 

17 

127.05 

1.  39 

9 

10.71 

3.00 

18 

85.54 

0.21 

TABLE  3.2 

Calculations  for  the  Optimum  Policy 


n 

xn+l 

yn+i 

Cn<TS» 

T* 

n+1 

cn«llTn.l’ 

2 

67.01 

1.56 

0 

46.73 

0 

37.72 

3 

65.47 

1.09 

0 

32.24 

0 

32.05 

4 

50.65 

0.91 

0.  33 

30.9 

0.  33 

27.06 

5 

47.85 

0.76 

0.92 

24.44 

0.92 

23.74 

6 

40.27 

0.69 

0.12 

23.00 

0.  13 

20.  55 

7 

34.17 

0.62 

0. 19 

21.43 

0.  19 

19.12 

8 

29.64 

0.77 

0.95 

25.6  3 

0.95 

23.  25 

9 

27.27 

1.02 

0 

26.21 

0 

24 . u \ 

10 

24.66 

0.93 

0 

26.23 

0 

24.32 

11 

27.27 

0.85 

0.09 

21.74 

0.  09 

23.63 

4.  CONCLUDING  REMARKS 


In  this  report  we  have  presented  two  models  for  the  determina- 
tion of  Bayesian  software  correction  limit  policies  under  the 
assumption  of  exponential  error  occurrence  times.  For  the  first 
model  we  assume  that  the  sampling  cost  is  negligible  while  for  the 
second  model,  such  cost  is  incorporated.  Procedures  for  determining 
the  optimum  policy  were  described  and  illustrated  via  numerical 
examples. 

I 


appendix  a 

TheoremA-1.  Let  C*(T)  be  a cost  function  given  by 

T 


C*  (T) 


a + C1  J G ( 1 1 £>dt  + c2  u2  G (Tl  £) 

t~3  I 

b + ; G (tl  £)dt  + n2  G(T  ly) 


(A-l) 


where  a and  b are  constants  and  let  the  following  conditions 
hold 


C1  < C2 


and  A 2 0 


where 


A = c^b  - a , 


Then  there  exists  a finite  and  unique  T*  which  satisfies 


q(T)  »r(TI^)|A+  (c2~cl)  J*  G(tl^)dt)  + (c  ,,-Cj  )G(TI  y)  = B 

where 


(A-2) 


B * — (c.b-a)  . 

^2  1 

Proof:  The  solution  of  the  equation  = 0 can  be  expressed 

in  terms  of  r(TI^)  and  is  given  by  equation  (A-2). 

Notice  that  the  repair  rate  r(TI^)  is  monotonely  decreasing 
with  T . It  can  be  easily  shown  that  the  IHS  of  the  above  equation, 
q (T ) is  also  monoton ically  decreasing  with  T under  the  conditions 
c1<c2  and  A*0  . Therefore,  if  q (0)  > R > q (<*>)  0 then  there 
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1 


exists  a unique  and  finite  solution  satisfying  the  above  equation. 
This  solution  yields  the  minimum  cost  because  the  cost  function 
is  convex  under  these  conditions.  Also,  from  the  monotonicity  of 
q(T)  , T*  * 0 if  B > q (0 ) and  T*  = ® if  B < q (®)  (see  Figure  A- 1 ) . 
If  there  exists  an  optimal  T*  then  by  substituting  T*  into  C*(T) 


we  get 


Cl-C2  ^2  r 

C*(T*  ) = l-4_  r (T*  ly) 


(A- 3) 


This  completes  the  proof. 


Note  that.  Theorems  2.1  and  2.3  are  the  special  cases  of  this 


general  theorem. 


r(Tly)  - (B-l) 

A+  (c2  - cj)  yn+i 

Proof  of  Theorem  B-l.  We  again  use  the  general  form  to  show  that 
there  exists  a unique  and  finite  upper  limit  T of  T*  in  case  of 
the  existence  of  T*  , i.e.  under  the  conditions  c^cc.,,  A»0 

and  q(0)  >B>q(®)  . Since  the  repair  rate  r(TI^)  is  monotonically 
decreasing  with  T , we  have 

r ( o * X ) > for  T>® 

or 

r(Ol^)  * G(TI^)  > g(TI^)  • 

Integrating  both  sides  over  the  range  of  T we  get 

*rH-lr(0,X)  > 1 (B-2 ) 

where 

*0+1  B J*0  G<T,j£>dt 

is  the  posterior  mean.  Now  define  the  following  function. 

M(T)  = r(TI^)[A+  (c2-c1)yn+1)  -q(T)  . (B-3) 

Then,  from  B-2  we  note  that  M (T)  is  monotonically  decreasing  in  T 
with 

M(»)  = 0 (B-4) 
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and 


M(0)  = <c2-Cl)yn+1 r(Oly) - (c2.Cl) 
= (c2-c1)(yn+1r(0ly)-l)>0 


Therefore, 


M(t)  > 0 for  T > 0 


or 


r(TI^)(A+  (c2~ci)yn+l^  >C*(T) 


Hence,  if  there  exists  a T*  such  that 


q (T*  ) = B 


(B-5 ) 


(B-6  ) 


(B-7  ) 


then  there  also  exists  a unique  and  finite  root  T satisfying 


r (T ly) { A + (c2_cl,yn+l^  = B 


or 


(T 


At  (Cj-Cjly,,,! 


(B-8) 


It  is  easily  seen  that  T > T*  (see  Figure  B-l) 


APPENDIX  C 

Relationship  between  C(T*)  and  T* 

n n 

If  c.  < c,  , then  C(T*)j  with  r(T*)  . This  is  seen  to  be  true 
from  equation  (2.25).  Now,  from  Figure  C-l,  we  see  that 

cnIT*>  » WT*>  <“>  > TS+l 

because  r(T)i  with  T.  Hence  the  decision  will  be  to  take  another 
observation. 

From  Figure  02,  we  see  that 

cn.i(T*> » v™  <_>  t;  •=  t;+i 

and  hence  the  decision  will  be  to  stop  taking  observations  and  to 
employ  the  optimal  policy  T*  . 
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